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Abstract 

The second-order partial derivatives of the Coulomb potential of a point charge can 
be regularized using the Coulomb potential of a charge of the oblate spheroidal shape 
that a moving rest-frame-spherical charge acquires by the Lorentz contraction. This 
'physical' regularization is shown to be fully equivalent to the standard delta-function 
identity involving these derivatives. 



Quantities with a singularity of the type at the origin r = occur in classical electrody- 
namics in connection with the idealization of a point charge distribution. For example, the 
straightforward calculation of the second-order partial derivatives of the Coulomb potential 
1/r of a unit point charge yields d^r'^ / dxidxj = {SxiXj — r^Sij)/r^, or the field of a point 
electric or magnetic dipole, obtained as the straightforward gradient of a potential with 
radial dependence has the radial dependence Because of the singularity at 

the origin, the integral of such a quantity over any three-dimensional region that includes 
the origin r = does not exist even in the improper-integral sense: the value of the integral 
obtained by excluding from the integration a region Vo around the origin and taking the 
limit of the size of Vq tending to zero depends on the shape and orientation of Vq. Integrals 
involving derivatives of or second-order derivatives of 1/r therefore have to be suitably 
regularized. A formal way of doing that is to use the delta-function identity [l| 



d"^ 1 3xiXj-r'^6ij 47r 

^k^ir) (1) 



dxidxj r 3 

where 6ij is the Kronecker delta symbol and 6{r) = 6{xi)6{x2)S{x3) is the three-dimensional 
delta function. The validity of the identity ((T)) can be justified most easily by the use of the 
straightforward regularization l/(r^ + a^)^/^ of the singular potential 1/r: 

lini ^ = lim ^^'^^ ~ ^^^'^ - lim (2) 

a^o dx^dxj + a2 a-o [r^ + a2)5/2 a^O (r2 + a2)5/2 ^ ^ 

since here the second term on the right-hand side is a well-known representation of 
— ^7i6ij6{r) (e.g., see 0). The first term on the right-hand side of the identity ((T)) is as 
such still non-integrable at the origin r = 0, but the regularization Q also includes a spec- 
ification of the regularization of this term; of course, the limits a ^ are understood to be 
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taken only after a three-dimensional integration with a well-behaved 'test' function. Regu- 
larizing the term (SxjXj— r^5jj)/r^ in a different, but equivalent, way, the identity may 
be written as 

9^ 1 3xiXj-r^6ij 4:71 
d^-r = - ^) - Y^--^(^) 

where 0(-) is the Heaviside step function. This formulation is equivalent to the stipulation 
that the spherical coordinates are to be used in the integration with a test function and that 
the angular integration is to be done first; the identity (0) has been derived in reference ji| 
effectively in the form Q. 

In a recent paper on the Coulomb-gauge vector potential of a uniformly moving point 
charge [3], an occasion has arisen of using the delta- function identity Q for the regularization 
of an integral of the type / d^r f{r)d'^r~^ /dxf in terms of which the difference between the 
Coulomb- and Lorenz-gauge vector potentials in that problem can be obtained as the solution 
to a Poisson equation. Since the relation (jS)) is an identity, there should be no doubt as 
to the correctness of such a formal regularization. However, in a problem that concerns a 
moving charge, it would be reassuring if one could show that a more 'physical' regularization 
procedure will yield the same results. Physically, it is natural to regularize the Coulomb 
potential 1/r of a moving point charge by the Coulomb potential of a charge that has 
the oblate spheroidal shape that a moving rest-frame-spherical 'elementary' charge of finite 
extension a acquires by the Lorentz contraction, and then to take the limit a — of any 
integral involving second-order derivatives of this potential. Such regularization involves an 
'ellipsoidal' approach toward the singularity at the origin because of the spheroidal shape 
of the moving elementary charge — in contrast to the first term on the right-hand side of 
identity © that stipulates a strictly 'spherical' approach toward this singularity. In this 
note, we demonstrate that a physical regularization along the above lines is indeed fully 
equivalent to the delta-function identity (jS)). We believe that the reward for carrying out 
the calculation that this demonstration requires will be a physical insight into and ensuing 
confidence in use of a useful formal relation. 

As a preliminary, we note that a 'physical' justification of the well-known delta-function 
identity V^(l/r) = —47i6{r) can be provided very simply. Let ^Pai'^) be the Coulomb poten- 
tial of a unit elementary charge described by the density Pa{r) = (1/Ki)0(a^— 7^x^— x^— 0:3), 
where Va = ^na^/'j, which is the density of a uniformly charged spheroid with semiaxes a/7, 
a, a, centred at the origin. For 7 = (1 — f^/c^)~^/^, this is the density of a charge that is 
moving with a velocity v (c is the speed of light) along the xi-axis and that is a uniform 
ball of radius a in its rest frame. Then, for any well-behaved test function /(r), we have by 
the fact that the potential (pa{f) satisfies the Poisson equation V'^ipa{r) = — 47rpa(r): 

/ dV/(r)VVa(r) = -An [ dV/(r)p.(r) = f dV/(r) = -|^K/(ro) (4) 

where the mean-value theorem is used on the right-hand side, with Tq being a point inside 
the region Va occupied by the spheroid. Taking now the limit a — 0, we obtain 

lim / dV/(r)VVa(r) = -47rlim/(ro) = -AnfU}) = -An [ d^rf(r)5(r) (5) 

because the point Tq G Va has to converge on the origin r = as a — > 0, and thus we can 
write 

limVVa(r) = -47r5(r). (6) 

a— »0 
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A uniformly charged ball was not the most popular model of an elementary charge em- 
ployed in the classical electron theory — this was a uniformly charged spherical shell, or, 
equivalently, a charged spherical conductor (see, e.g., jj]). The Coulomb potential of a uni- 
formly moving charged conductor that is spherical with radius a in its rest frame is the most 
convenient one to use for our purpose because it equals the Coulomb potential of a charged 
conducting oblate spheroid of semiaxes a/7, a, a (see j^f; an interesting historical back- 
ground to the problem of a moving charged sphere can be found in j^), and this potential 
can be expressed in terms of an elementary function Fzl] : 



{l/(3a) arctan[/3a/^i(r++r_)2-(/5a)2] for -i^xl + > a"^ 
{1/ 13a) arctan(/37) for 7^0;^ + < a? 



where 



r± = ^xl + {p±(3aY p=\Jxl + xl 7 = -^=L=^ /5 = ^. (8) 



The partial derivatives d^pair) j dxi of the potential '^aix) for ^ x\ + p >a are 



dxi 



Xl 



dx 



2,3 



{r_ -r^)l3a+ {r^ + r_)p 0:2,3 



(9) 
(10) 



Since the potential f^a^f) is constant inside the spheroid '-y'^x\+p^ = a?, the partial derivatives 
dipa{r) / dxi can be written as Q[{'j'^xl+p'^Y^'^—a]d(faif)/dxi, and thus the second-order 
partial derivatives d'^ipa{r) / dxidxj can be written as 



n n . n e(v/7%?+P^-a) + -^^Q{^l'xl+p^-a) 

OXiOXj OXiOXj ^ OXi OXj ^ 

9^'^a{r) ^, / 2 2 , 2 ^ , ^V^air) [l+{-f^-l)6ij]xj / , 2 , 2 

+ ~d^- / o , . o ' '^(\/7'3;f+p2-a) (11) 



^72x? + p2 



where the derivatives of the potential on the right-hand side are understood as those of the 
expression for the potential exterior to the spheroid. 

The equivalence of the delta-function identity Q and the regularization that uses the 
Coulomb potential (pair) of a charged conducting spheroid demands that, for any well- 
behaved test function /(r). 



d^^a{r) f .3 3xiXj-r^6ij An 
mil I , J )— — - — = lim / d r/(r) ^— ■ 

a^Oj OXiOXj 0.^0+ Jr>a 3 



hm^ dV/(r)-^ = hm dV/(r) "-^-^^ ^ - ^5.,/(0). (12) 



Using (fTT|) . we write the left-hand side of (fT^ as 

lim /dV/(r)^^=lim / dV/(r)^^ 

a-^Oj OXiOXj a-*0Ji^xl+p'^>a2 OXiOXj 



2 iu ir /^s , gv^a(r-) x,-(5(v^72x2+p2-a) 
+ [l+(7 -l)5i,]limy d r/(r)^- ^_=^.(13) 
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Using the expressions ^ and (jTUl) for the derivatives dipaif)/dxi, the second term on the 
right-hand side of ()13|) can be evaluated in cfosed form. Let us first assume that i=j=l. 
Transforming as 7x1 Xi and then to the spherical coordinates r, 6, 0, with xi as the polar 
axis and cos^ = ^, we obtain 



7^ lim / d^r f(r) 



^-1^X1 + p2 



= lini / dV/(xi/7,X2,X3; 



Xi5(r — a) 



27r ^2 



/ d0/(Xi/7,X2,X3)|^=a 
-1 JO 



./ ^ f2 / 2 2arcsin/3\ , 

where dipa{r) / dxi\xi^xi_h denotes the exterior partial derivative © after the transformation 
7x1 — > Xi- Here, the delta function 5{r — a) led to an immediate radial integration, which 
enabled a considerable simplification of the integrand; the limit a — then could be taken 
inside the remaining integral, yielding 

lim /(xi/7, X2, X3)|r=a = lim /(a COS 6^/7, a COS sin 6*, a sin sin 6*) 

a— >0 a— >0 

= /(0, 0,0)^/(0). (15) 
A similar calculation for i=j=2 yields 



lim / d^r /(r' 
a— >o J 



difair) X25(^7^x? + p2_a) 



9X2 ^7^2 + p2 

X2S{r — a) 



lini / /(a;i/7,X2,X3^ 



a-^O J 7 dX2 

/(o) ^^^^ ^Ai-e)cos'^ 



72 7-1 ^1-/32(1-^2) 

-r = -27r 1- — + /(O). (16) 



72 7-1 M - /?2(1 _ ^2) \^ ^2 ^^3 

Here, similarly as in the integration in (fT^ . even with the relatively complicated expression 
(fTUI) for the exterior partial derivative dipa{r)/dx2, the integrand could be simplified con- 
siderably after the radial integration. The case i=j=3 will obviously yield the same result, 
while for any mixed case i j, the integration with respect to the azimuthal angle (p will 
lead to a vanishing result. Collecting these results, we have 



a^UJ aXi ^^2j^2_^p2 

where 

2//?2 - (2/7/53) arcsin/3 for i=j=l 
^/,,(/3) = <! 1 - 1//?' + (1/7/?') arcsin/5 for t=j=2, 3 (18) 
for i 7^ j. 
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We note that lim^^o = ^^ij- Using (fT^ and (fTTj) . the condition (fT^ of regularization 

equivalence can now be written as 

li- . _ d\ fir) - hm / dV /(r) ^^^^^^ " 



^5 



a^0 7722;2+p2>a2 OXidXj 0.^0+ Jr>a T 

= 2n[g,,i^) - f<5,,-]/(0). (19) 

To prove that the condition (fT^ holds true, we proceed as follows. The first limit on the 
left-hand side of ()19|) can be written more simply as 



lim / dV/(r) = lim / d'r f(r) ' ' , — (20) 

1^0 7723,2 +p2>a2 dXiOXj a->0 7^2^.2+^2 >a2 

since the derivatives d^ipa{r)/dxidxj for 7^x^ + > can be expanded in powers of (/3a)^ 
(such an expansion is obtained most easily by differentiating term by term the corresponding 
expansion of ipa), where the zeroth-order term equals [SxiXj — r'^6ij)/r^ and the integrals 
involving the higher-order terms vanish in the limit a —>■ 0. It suffices to show this only for 
the case i=j=l. Here we have 



5Va(r') 



dxi 



^ , = ti-m^n + 2)P,„,.(x,/r)^^, (21) 

72x^+p2>a2 I 



where Pm{-) are the Legendre polynomials. Using (PT|) and the multipole expansion of the 
function f{r) after the transformation Xi —>■ Xi/j, 

/(r cos 6^/7, r cos sin 61, r sin sin 6^) = ^ /^^(r, 7)1^/^(6', 0) (22) 

Im 

we obtain 

, ,dV/(r)%M = f e„(a,7)(/5a)^" (23) 

where 

//m(r,7) 

Im 

P2n+2icOs9/'Ju) 



c^{a, 7) = E CnUl) r dr (24) 
Cnimil) = (-l)"(2n + 2) J d^]l^„(^,0)^^^^^±^^^^^/^ u = ^l-(3'cos'e. (25) 



We note that C„oo(7) = for any n > 0. As limr-+o /«m('", 7) = for any / > 0, we 

have that lima^ofa^" XT' dr /^^(r, 7)/r^"+^] = when both n > and I > 0, and thus 
lima^o [cn (a, 7) (/^o)^"] = for any n > 0. Therefore, we indeed obtain 

lim/ dV/(r)^^^4^ = limco(a,7) = lim / dV/(r)^^i^. (26) 

a-*0y72^2+p2>„2 ' ^ ^ dxj a->0 ' a^O ^2^2_^.p2>„2 ^ ^ ^5 ^ ^ 

Using (pUj) . the condition (fT^ can be expressed as 

lim / dV /(r) ^^-^^ :"'^-^- = 2n[g^,iP) - f 5.,]/(0) (27) 



6 



where the integration region Ua is the region between the surfaces of the oblate spheroid 
7^x^ + = a? and the sphere x\ + = a^: 

Ua = {(xi, X2, Xs); 7^x^ + x\ + x\> a^, x\ + x\ + x\ < a^}. (28) 

We now evaluate the left-hand side of (j2Z|)- When the size parameter a tends to zero, the 
integration region Ua gets progressively smaller and closer to the origin r = 0, and thus, for 
r G Ua, f{r) — > /(O) as a 0. We can therefore write the left-hand side of as 



lim / dV/(r)^^^^^-^ = /(0)lim / dV^^^^^^^. (29) 

Transforming here the integral on the right-hand side to the spherical coordinates, with xi 
as the polar axis and cos^ = ^, we obtain for i=j=l: 



|'^de(3e^-l)ln[l + (f -l)e^] 



„ /2 2arcsin/? 2\ , . 

2^ - - o • (30) 



/32 7/33 3 



The case i=j=2 gives 

- 1 /•! _ /•27r 



J lAa ^ Zi J — ^ 1 >/ 

£d^(i-3aMi+(f -m'] 



TT 

2 



/I 1 arcsin/5\ 

= '''' 

and the same result will obviously be obtained for i=j=3. The mixed cases i ^ j will 
all yield zero on account of the integration with respect to 0. The values of the integrals 
(jHU)) and are independent of a, and using these results and (|^. we obtain (j2Zj). This 
completes the proof of the regularization equivalence (fT^ . 

In closing, we would like to stress that the regularization equivalence p2p is bound to 
hold also when the potential ^Pa{f) is the Coulomb potential of a uniformly charged spheroid, 
or of any other Lorentz-contracted charge distribution that is, in its rest frame, spherically 
symmetric and characterized by a finite size parameter a. (Explicit expressions for the 
potential of a uniformly charged spheroid can be found in the literature 0, |^ [l^ , but they 
are rather more complicated than the conducting-spheroid expression ((7j).) In fact, using the 
powerful results of the theory of generalized functions and derivatives, one can very easily 
give a formal proof that ()12j) holds for the Coulomb potential fai^) of a charge distribution 
Pa{r) of any shape, subject only to the condition lima_>o Pa(^) = this proof is given in 
Appendix. 
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Appendix 

Lemma. Let 



^3 / Pair') 



^«W=ydV^, (32) 

where the integration extends over all space and Pa{r) is a localized (not necessarily spher- 
ically symmetric) function of r = {xi, X2, x^) that depends on a parameter a so that 

wlim pa(r) ^ 5 (r). (33) 

a— >U 

Then 

whm = -—5,,Sir) + w lim ^ 9 r - £ , 34 

a^o dxiOXj 3 £^0+ 

where r = |r| = (x^ + 2^2 + and ©(■) is the Heaviside step function. The symbol wlim 

denotes the weak limit: wMuia^ao fair) = g{r) iff lima^a^ / d^rt{r)fa{r) = J d^rt{r)g{r) 
for any 'well-behaved' test function t{r). 

Proof. To prove (34), we need to show that, for any 'well-behaved' test function /(r), 

lim / dV f{r)^^ - -^fmj + lim / dV f{r) ^^^^^ "'^'^ (35) 

To evaluate the left-hand side of (35), we note that we can replace the derivative d"^ /dxidxj 
by the generahzed (distributional) derivative B^/dxidxj. This allows us to exchange the 
order of the hmit a and the differentiation [11, p 12] since the space of the generalized 
functions is complete [12]: 



a-^Oj •'^ ' dXidXj a^oj ^ ' dXidXj 



92 



j d. r f{r) Q^^Q^^ ^ ■ (36) 



Here, the 3rd line was obtained using (32) and (33). But 



1 47r ^ ^, , 2)XiXj — r'^SiT ^ , , 
— — = -—5i.5ir) + w lim ^ , ^ 6 r - £ 37 

(see [11], p 28, but note that the signs of the right-hand sides of (3.129) and (3.130) are 
there misprinted; see also [13,14]), and thus 

Using (38) in (36) results in (35). QED 

The author is grateful to F Farassat for a useful discussion on the matter of this appendix. 
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